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Abstract: We consider large N gauge theory on a two dimensional lattice in the presence 
of a baryonic chemical potential. We work with one copy of naive fermion and argue 
that reduction holds even in the presence of a chemical potential. Analytical arguments 
supported by numerical studies show that there is no phase transition as a function of the 
baryonic chemical potential. 
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1. Introduction 

The ! t Hooft model jy], namely, large iV gauge theory with finite number of fermion fla- 
vors in two dimensions, has several similarities with four dimensional QCD. The particle 
spectrum has an infinite tower of mesons [Q], and there is a massless meson in the chiral 
limit. Chiral symmetry is broken in the N —* 00 limit even though it is a two dimensional 
model. Confinement is simply realized since Wilson loops obey an exact area law [Q, Q]. 

Eguchi-Kawai reduction |3j holds in two dimensions || and the model can be reduced 
to a single site on the lattice. The two Zn symmetries (U(l) in the N — * 00 limit) 
associated with the Polyakov loops in the two directions remain unbroken as one approaches 
the continuum limit. Therefore, the model does not depend on the spatial extent or the 
temperature and it is always in the confined phase. Since finite number of fermion flavors 
do not affect the gauge field dynamics in the N — > 00 limit ||, the chiral condensate can 
be computed in the "quenched approximation" and is independent of the temperature. 

We will study the role of the baryonic chemical potential in the Euclidean formalism of 
lattice large N QCD. The Zn symmetries associated with Polyakov loops will play a central 
role in our discussion and we will use a generalized partition function to obtain a result 
for the quark number density with a single copy (four degenerate flavors) of naive fermion. 
The main result of this paper is the independence of the physics on the quark chemical 
potential. We provide analytical arguments and support it with numerical calculations. 
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2. The quark number density 



We consider the generalized partition function of lattice QCD on a two dimensional L\ x L 2 
lattice with one copy of naive fermion defined by 

z(4>,x) = [ n[^i^)][^2(x)]e 2W S,Re^, det ^(r 1 ,r 2 ,m ff ,^x), (2.1) 

where U fJi (x) are SU(N) matrices on the links connecting site x to x + ft for fi = 1,2. 
[<ZZ7i(x)] and \dU2{x)\ are the standard Haar measures on SU(N). b is the lattice gauge 
coupling which is kept constant as N is taken to infinity. The continuum limit is obtained 
by taking b — ► 00 and the lattice spacing can be set to U p is the standard plaquette 
operator. <p and x are two independent complex variables. We will assume that iV is odd 
so that the baryons are made up of an odd number of quarks. The naive Dirac operator is 

M {T Xl T 2 ,m q , cf>, X ) = 2m q + a x - ^T/) + a 2 ( X T 2 - x" 1 ^) , (2.2) 

where are unitary operators given by 

[T^]{x) = U PL {x)^{x + ji). (2.3) 

M is a matrix with a linear extent of M = 2NV with V = L\L 2 - 
It is clear that the result is a finite polynomial of the form 

M M-\k 2 \ 

Z (^X)= £ £ a klM ^ X k2 - (2.4) 

k 2 =-M ki=-M+\k 2 \ 

Since the result is a finite polynomial, one can obtain all the coefficients by setting eft = e ip 
and x = elul with p, uj G [— tt,tt} and using fourier transforms. This is by no means a new 
idea || nor does it avoid the sign problem since it involves fourier transforms of measured 
quantities. 

The fermion determinant appearing in the definition of Z(e ip ,e lul ) is positive definite 
and it does not affect the gauge dynamics since the effect of p and uj is simply to make 
the fermions interact with a constant external U(l) field on top of the SU(N) gauge field. 
Therefore, we can invoke Eguchi-Kawai 



X 

jfdetX (^ 2 ,m g ,e< p+ ^),e<" + ^r)) . (2.5) 

— n l n >. / 



fc 1= k 2 =0 

We expect factorization of observables in the N — > 00 limit and therefore, we can write 

Zf(e ip e tuJ ) = { ' 

/nj^i][^2]e 2 ™ eWlC/s[/ '^ 



1 This amounts to fixing the string tension at 4 as one approaches the continuum limit 
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Lj-lLa-l 



II II 2?/ ( e^ p+ ^fi , e^+T^ ) (2.6) 

fe 1= k 2 =0 



where 



v ) = I n,[^i] [dU2}e 2bNRe ™^fM det Aj (Eft, Eft, m,, 0, x ) 



The gauge action and the Haar measure in (E.7[) are invariant under C/i — > e at £/i and 

■2vr — 

[/ 2 — ► e n JJ 2 . Therefore, 

Vf{tf>,x)=V f {f^<f>,xj =T> f (<l>,e i %x) V tf>, X . (2.8) 



It follows from (fO|) that 

o-i^(f/i,E/ 2 ,m ? ,^,x)o-i =M(U 1 ,U 2 ,m q ,(f),-x) (2.9) 

and 

a 2 ^4(C/ 1) [/ 2 ,m 3 ,^x)^ 2 = >t(^i,^2,m 9 ,-0 ! x). (2.10) 

Equations (E^§1), (l2.9|) and ( 2.1C| ) along with (0) restricted to V = 1 results in a 
polynomial of the form 

Vfifa x) = a ,o + a 2 N,o<f> 2N + a_2Ar,o0 _27V + a ,27VX 27V + a -2NX~ 2N ■ (2-11) 

Since a±2jv,o and ao,±27V correspond to the highest or lowest term in the polynomial in 
the respective variables, it follows that 

a 2 N,o = (detaiUi); a- 2 N,o = (det o\U{); a , 2 N = (det a 2 U 2 ); a = (det a 2 Ul), 

(2-12) 

and hence all these coefficients are —1 for odd N. 
Therefore, we finally have, 

2>/(6 x) = ao.o - (4? N + r 2N ) ~ (X 2N + X~ 2N ) (2.13) 

The number density, p(fi), at zero temperature and infinite spatial extent (L±, L 2 — > oo) 
at a fixed chemical potential, /i, is given by 

/>(m) = I™ — ^— |~logZ(l,^) (2.14) 
Li,i 2 -»cx3 NL\L 2 op 



Using (P^q) and ( |2.13|) , we get 

1 /' 7r 

p(p) = —j dpn(p,p) (2.15) 

where 

I , -Usm.2N(uj - ip) 



n(p, p) = — / da; ; 7. (2.16) 

KF,r> 2ir J_ n a 0fi - 2cos2Np - 2cos2N(co -ip) y } 
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This result for the number density with an explicit integral over the momenta, p, and the 
Matsubara frequencies, u, for the interacting theory in the large N limit is an example of 
the quenched momentum prescription ||. 
Let us set 



a , = 4M 2 + 4. 



We can rewrite (2.16) as 



n{p, 11) 



AN Fn 



ITT 



sin 2N (u — ifi) 



AM 2 + 2 + 4 sin 2 Np - 2 cos 2N(uj - z» ' 



(2.17) 



(2.18) 



With z = e 2NluJ + 2N v ; th e above integral becomes a contour integral over a circle of radius 
e 2Nfi cen t erec i a t the origin and 



n(p,fi) = — 

ITT 



1 



z [z 2 + 1 - 2 (1 + 2 {M 2 + sin 2 iV>}) 
if n < sintT 1 yW 2 iVp + M 2 
2 if /x > i sinh- 1 ^sin 2 Np + M 2 



(2.19) 



Writing, fj, = J±= with /i being the physical chemical potential, we see that the critical 
chemical potential, p, c , is given by 



u c = lim \fb lim — sinh 1 

b^oo N— >oo TV 



sin 2 Np + M 2 . 



(2.20) 



By comparing this result with that for free quarks, namely, 



»KP>A*J in f z [ z 2 + i -2(1 + 2 {m 2 + sm 2 p}) 
if fi < sinh -1 * / sin 2 p + m 2 



-dz 



2 if n > sinh 1 \j sin 2 p + m 2 



(2.21) 



we see that the effect of the interactions is to scale \x and p by a factor of N (Np and iV/i 
can be viewed as the baryonic mometum and baryonic chemical potential) and replace the 
quark mass, m q , by M defined as 



M 2 = jD f (l, 1) = i (det 2m g + a x [ U x - U{ ) + <r 2 ( *7 2 - U. 



(2.22) 



3. The physics of M 



It is clear from ( 2.22 ) that M = m q for f/i = U2 = 1- It is also clear that M will behave 
like e QiV for large iV and a will depend on b and m 3 . Following 't Hooft we will keep 



7 = 2irbm 2 fixed as we take the continuum limit, b 



00. 



We can compute the "tree level" contribution to M by setting 



e idiag(«i,« 2) — >ujv) jj 



idiag(«i,i;2,— ,«jv) 



(3.1) 
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where Ui,Vi G [— vr,7r] and u i = Yld=i v i = since we are in the confined phase. 

Setting m q = 0, we have 

M 2 = Jjv(0,0), (3.2) 
where In(p, q) is defined in Appendix [A| with 

f(u, v) = A (sin 2 u + sin 2 v) = 4 - (e 2iu + e" 2 ™ + e 2 ™ + e" 2 ™) . (3.3) 

Therefore, 

M 2 = A N + A{-1) N , (3.4) 

and the critical chemical potential p, c in ( |2.20| ) will approach infinity as ln(2)v6 when 
b — > oo. This is the main result of our paper. 



Following the steps similar to the ones from (2.14) to ( 2.19| ), we can show that the 
chiral condensate, x(m)> defined as 

XM= lim lim — — 3 logZ(l,e") (3.5) 

m q ->0Li,L 2 ->oo N L\Li om q 

reduces to 



[0 if \i > ^sinh" 1 ^/ sin 2 Np + M 2 

Since the chiral condensate should behave as 



X(0) = (3.7) 



we expect a to take the form 



as we approach the continuum limit for large N. 
4. Numerical estimate of a 

We support our argument in Section |3| by performing a numerical estimate of a. The main 
purpose is to convincingly show the presence of ln(2) in ( |3.8| ). Our numerical procedure 
has the following steps: 



1. We fix N, b and 7 and numerically evaluate the right hand side of ( p. 22 ). 



2. We plot |ln2?j(l,l) — iVm2 as a function of N at a fixed b and N to extract the 
large N limit of (a — ln(2)). 



3. Based on (|3.8|) , we fit 

a-ln(2) = Co(6) + C 1 (6)^7 + 0( 7 ) (4.1) 
and extract Co (6) and C%(b). 
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4. We fit C (b) to 

C = zo + ^ b + r f+0(b-^) (4.2) 
with the aim of showing that zq is consistent with zero. 

5. We fit d{b) to 

with the aim of showing that z\ and [i\ are consistent with zero. 

Our primary aim, as mentioned above, is to show that Zo, z\ and [i\ are consistent 
with zero. The gauge action was simulated using a Hybrid Monte Carlo algorithm since it 
performs better than a Metropolis algorithm based on SU(2) subgroups. The simulations 
were performed for all odd values of N in the range {11, 39}. For each N we used 12 values 
of b = 3, 4, 5, 6, 7, 8, 9, 14, 20, 40, 60, 100. We focused on going to smaller values of b in order 
to get close to the tree level result in Section [| The HMC time step ranged from 0.02 
(at small b) to 0.004 (at large b). The number of steps in a trajectory was chosen so that 
the length of the trajectory was 1.0. The fermion determinant on each configuration was 
obtained by performing an average over N 2 transformed gauge field configurations: 
e l ^~Ui,e l ~N~U2)] < k\,ki < N. This numerically enforces (p^q). We computed 
the fermion determinant for 21 different values of masses in the range 7 = [0, 4] on every 
configuration. 

A sample plot at b = 3, 8, 20, 100 and 7 = 3 shown in Fig. || illustrates the extraction 
of (a — ln(2)). It is clear from the slopes in Fig. |l| that m(2) dominates the value of a 
in our range of b. Yet, we have sufficient statistics to extract (a — m(2)) with reasonable 
accuracy even at the weakest coupling. 

A sample plot at b = 3,8,20, 100 shown in Fig. Q illustrates the extraction of Co(fe) 
and Ci(b). The shaded data points with 7 > 1 were used to perform the fit. Deviation 
from the linear behavior for 7 < 1 is a consequence of finite N effects. 7 = 1 is roughly 
the transition point between light and heavy quarks |lC| l and therefore we do not expect 
to obtain a good estimate of S. One should note that the errors in the intercept, Co (6), 
are larger than the errors in the slope, C\(b), and this is a consequence of subtracting the 
dominant ln(2) term. 

The plot of Co (6) vs ^= is shown in Fig. ||. The first observation is that the whole 
range of the plot is less than 0.21n(2). The second observation is a clear sign that Co(b) 
tends to zero as b — > 00. The fit of Co (b) to ( f4.2p , shows that zq is zero within the expected 
accuracy of this quantity. The other two fits, one with zq forced to zero and the other with 
Ho forced to zero, show a larger fluctuation in the non-zero coefficients indicating that we 
do not have enough data to extract /io or r^. 

The plot of Ci(6) vs ^= is shown in Fig. |j. A fit to the form in ( |4.3| ) is shown with 
solid lines. The fit is consistent with z\ = and [i\ = 0. To further emphasize this 
point, a fit with z\ forced to zero shows that [X\ is still consistent with zero. Finally, a 
comparison of these two fits to a fit with only the \ term shows that the coefficient of the 
non-zero term remains unchanged in the three different cases. Since naive fermions in two 
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dimensions describe four flavors of fermions, the condensate for a single flavor is estimated 
as | = 0.34(3). This estimate is significantly higher than the analytical result [^] but this 
is to be expected since we only used 7 > 1 in our fits. 



5. Discussion 

The analysis in section || shows that the absence of a critical chemical potential in two di- 
mensional large N gauge theory can be seen at the tree level. Numerical analysis presented 
in section || supports the tree level argument. The basic result is that the critical chemical 
potential behaves as (ln2)\/fr for naive fermions and one can extract the coefficient of the 
term using a tree level calculation. The coefficient of the y/b term (equivalently, the 
constant term in (|3.8| )) will depend on the fermion discretization but one can show on 
general grounds that this coefficient is positive for any discretization. The doubly periodic 
function, f(u,v), in ( |3,3[ ) is the naive fermion determinant on a single site lattice in the 
presence of a U(l) gauge field of the form (e tu , e w ) and is a positive function for all (u, v). 
Different fermion discretizations will lead to different expressions for f(u,v) but they will 
all have the general property that it is a positive function for all (u,v). It follows that 

7(0, 0) > 7(n, m) V(n, m) + (0, 0) (5.1) 
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in ( |A.1| ). Since /(n, m) are raised to the A^th power upon integration over the gauge 
fields in ( |3.1| ) as a consequence of the identity in Appendix [A], it follows that 7^(0,0) will 
dominate in the large A" limit resulting in a positive constant term in ( |3.8| ) and there by 
driving the critical chemical potential to infinity. 

The above argument shows that there is no critical chemical potential for quarks in 
two dimensional large A' QCD. If Eguchi-Kawai reduction is not valid, as is the case in 
d > 2, then we cannot work on a finite lattice in the continuum limit and the arguments 
in this paper do not hold. 



Two dimensional QCD has also been studied in the Hamiltonian formalism [11, 12] 
Baryons have been discussed in the Hartree-Fock approximation |ll[ and the role of chem 



ical potential and temperature have also been discussed [12|. Contrary to the Euclidean 
Lagrangian formalism where reduction makes the independence on temperature quite ob- 



vious, it is harder to see the same result in the Hamiltonian formalism [12|. Arguments 
in support of the existence of a critical chemical potential is provided in [^] along with 
the possible breakdown of translational invariance at high densities. The dependence on 
spatial momentum is suppressed in our analysis since M dominates the denominator of 
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Figure 3: Plot of Cq(6) as a function of b. 



the integrand in Q2.18|) . This can also be seen in the expression for the critical chemical 
potential in ( 2.20|) and this is a consequence of the integration over the Polyakov loop vari- 
ables defined in ( |3.l[) . The Hamiltonian formalism, by construction, is in the continuum 
time limit and starts with the theory in the Weyl (temporal) gauge (Aq = 0) and imposes 
periodic boundary conditions on the spatial gauge field in the time direction. This amounts 
to setting the Polyakov loop in the time direction to unity which is not consistent with the 
unbroken Zfq symmetry. Recently, the importance of the role of Polyakov loop variables 
in the spatial direction has been addressed in the Hamiltonian formalism [in 
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A. A useful integral 

Let f(u,v) be a periodic function of u and v with period 2tt. We decompose it into its 
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fourier components according to 

f( U ,v)= y, E /( 

r 

Then, the integral 



oo oo 



n=— oo m=— oo 



n, m le 



t(reu+m») 



r N , N , N / N \ / N 

= / n ? n P n / («* + * + 9) 2 ^ e £ 

^ i=i 71 i=i ^ i=i \i=i / \i=i 



is 



00 00 



%0,g) = E E 



n;=— 00 mi=— 00 
AT 



JV 



/(n is 77^) 



i=i 
A 7 " 



(A.l) 



(A.2) 



J i=i J i=i \i=i / \i=i 



i=i " i=i 
00 00 

£ £ f N (n,m)e 



iN(np+mq) 



(A.3) 



n=— 00 m=— oo 



- 10 - 



References 

[1] G. 't Hooft, Nucl. Phys. B 75, 461 (1974). 

[2] A. A. Migdal, Sov. Phys. JETP 42, 413 (1975) [Zh. Eksp. Tcor. Fiz. 69, 810 (1975)]. 

[3] B. E. Rusakov, Mod. Phys. Lett. A 5, 693 (1990). 

[4] T. Eguchi and H. Kawai, Phys. Rev. Lett. 48, 1063 (1982). 

[5] G. Bhanot, U. M. Heller and H. Neubcrgcr, Phys. Lett. B 113, 47 (1982). 

[6] G. 't Hooft, Nucl. Phys. B 72, 461 (1974). 

[7] D. J. Gross and E. Witten, Phys. Rev. D 21, 446 (1980). 

[8] M. P. Lombardo, PoS C POD2006, 003 (2006) [arXiv:hep-lat/0612017]. 

[9] D. J. Gross and Y. Kitazawa, Nucl. Phys. B 206, 440 (1982). 

[10] J. Kiskis, R. Narayanan and H. Neubcrgcr, Phys. Rev. D 66, 025019 (2002) 
[arXiv:hcp-lat /0203005] . 

[11] L. L. Salcedo, S. Levit and J. W. Negele, Nucl. Phys. B 361, 585 (1991). 

[12] V. Schon and M. Thics, arXiv:hep-th/0008175. 

[13] B. Bringoltz, arXiv:081 1.4141 [hcp-lat]. 



- 11 - 



